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We assess the importance of final state interactions in — >■ K~ii^-k'^ , stressing the consistency 
between two- and three-body interactions. The basic building block in the calculation is a Kn 
amplitude based on unitarized chiral perturbation theory and with parameters determined by a fit 
to elastic LASS data. Its analytic extension to the second sheet allows the determination of two 
poles, associated with the k, and the ii'*(1430), and a representation of the amplitude based on 
them is constructed. The problem of unitarity in the three-body system is formulated in terms of 
an integral equation, inspired in the Faddeev formalism, which implements a convolution between 
the weak vertex and the final state hadronic interaction. Three different topologies are considered 
for the former and, subsequently, the decay amplitude is expressed as a perturbation series. Each 
term in this series is systematically related to the previous one and a re-summation was performed. 
Remaining effects owing to single and double rescattering processes were then added and results 
compared to FOCUS data. We found that proper three-body effects are important at threshold 
and fade away rapidly at higher energies. Our model, based on a vector weak vertex, can describe 
qualitative features of the modulus of the decay amplitude and agrees well with its phase in the 
elastic region. 

PACS numbers: 13.25.Ft ll.80.Jy 13.75.Lb 



I. MOTIVATION 

About forty years ago, reactions of the type KN — > 
ttKN were used to determine the Ktt amplitude [1] . Such 
reactions involve the scattering of an incoming kaon and 
a pion from the nucleon cloud. The dominant one-pion 
exchange ampUtude is isolated by selecting events with 
low momentum transfer. This is the only Ktt — > Kn 
scattering data, collected in the range 0.825 < mKir < 
1.960 GeV/c2. 

In the past decade, heavy flavor decays, in particu- 
lar decays of D mesons, became a key to the physics of 
the light scalars. Currently these are the only process 
in which S-wave amplitudes can be continuously stud- 
ied, starting from threshold, filling the existing gaps on 
the scattering data. In addition, very large, high pu- 
rity samples, in which the initial state has always well 
defined quantum numbers, became available in the past 
few years. Multibody decays of heavy flavor particles 
proceed almost entirely via intermediate states involving 
resonances that couple to tttt and Kn. The universal 
Kit and tttt amplitudes are, therefore, present in these 
decays as well. These amplitudes could, in principle, be 
extracted with great precision. 



* patricia@if.usp.br 
t tobias@ita.br 
alberto@cbpf.br 



Most of the existing results come from hadronic de- 
cays of D mesons. The golden modes are the D^,Df — >■ 
TT~7T~^7T~^[2] and, in the case of the Ktt system, the 
— > K^TT^TT~^[3] decay. These decay modes share 
some common features: the presence of two identical par- 
ticles in the final state, and a largely dominant S-wave 
component. The standard procedure is the analysis of 
the Dalitz plot, in which the decay amplitude is modeled 
by a coherent sum of resonant amplitudes, accounting 
for the possible intermediate states - the so-called iso- 
bar model. The extraction of the resonance parameters 
and decay fractions, however, depends strongly on the 
particular model used for the S-wave. 

The situation concerning the experimental results on 
the Ktt amplitude is intriguing. The S-wave amplitude 
can be also studied with semi-leptonic decays. In prin- 
ciple, in decays such as D ^ KttIvi or t — ;> KttVt, 
the extraction of the S-wave would be simpler than in 
the case of hadronic decays, for the Ktt system is free 
from final state interactions (FSIs) with the lepton pair. 
The Ktt S-wave, in this case, should match that of 
LASS[1], provided no energy- dependent phase is inher- 
ited from the weak decay vertex. The results are con- 
flicting, though. While the BaBar analysis of the decay 
D+ K-TT+e+Ve [4] and the FOCUS analysis of the 
decay — > K^tt^ fi^i/^ [5] conclude that the Ktt am- 
plitude is consistent with the LASS results, the analysis 
of the decay — > K'^tt~i't, carried out by the BaBar 
and Belle Collaborations [6, 7] showed that these data 
cannot be described by the LASS amplitude. 
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In 2006 the E791 Collaboration published a model in- 
dependent analysis of the K~'k'^ S-wave amplitude using 
the K~Tr~^TT^ decay [8]. A very similar analysis 

was performed by the FOCUS Collaboration [9]. The 
CLEO-c Collaboration also studied this decay, but with 
a somewhat different method [10]. In the E791 approach, 
the S-wave K~n~^ amplitude is represented by an un- 
known complex function of the K~Tr~^ mass, to be deter- 
mined directly from the data. The P- and D-wavc were 
parameterized by the usual sum of Breit-Wigner ampli- 
tudes. The /r~7r+ mass spectrum was uniformly divided 
into 40 bins. In each bin the S-wave amplitude was de- 
fined by two real numbers, Ao{m-'j(^) = aQe^^\ The set 
of 40 pairs (oq, S^) (80 free parameters) define the S-wave 
through the entire Kn spectrum. The phase and magni- 
tude of the S-wave at an arbitrary position were obtained 
by a cubic spline interpolation. 

The K~TT'^ S-wave amplitude obtained by E791 and 
FOCUS is significantly different from that from LASS. 
Possible explanations for this discrepancy fall into two 
broad categories and the first one concerns the weak ver- 
tex. In the hadronic scale, the mass of the D is far both 
from chiral and heavy quark limits, and the treatment 
of its decays cannot be simplified by the use techniques 
developed in these realms [11]. The second class of ef- 
fects concerns strong interactions, which do take place 
after the weak decay. The treatment of this part of the 
problem is necessarily involved since, even in its simplest 
version, these FSIs already involve three bodies (sec e.g. 
[12]). Other examples of the importance of final state in- 
teractions in three body decays of heavy mesons can be 
found in [13-15] As the final mesons are light and have 
high energies, kinematics is fully relativistic and tech- 
niques developed in low-energy Niiclear Physics, for the 
treatment of three-body systems, do not apply. In the 
case of relativistic Particle Physics, in spite of a grow- 
ing literature [16, 17], the corresponding techniques are 
still in the want of being developed for the application to 
the decay problem. Approaches to the relativistic few- 
body systems (two and three-body) have been collected 
in a series of works that presents the main developments 
done so far (see references [18-24]). The weak vertex of 
D+ ~> if ^7r+7r+ was treated by means of factorization 
and form factors, supplemented by two-body final state 
interactions [25, 26]. 

In this work we concentrate on the three-body struc- 
ture of strong dynamics of FSIs and postpone a detailed 
discussion of the weak vertex. As one knows little about 
this problem, our aim is to identify leading efl:ccts and 
a number of simplifications are made. Hopefully, once 
these leading effects are understood, corrections can be 
included in a systematic way. 

Our discussion of three-body FSIs is based on the the- 
oretical framework provided by chiral symmetry in the 
SU{3) sector to describe the two- meson scattering ampli- 
tude. At low energies, chiral perturbation theory (ChPT) 
yields the most reliable representation of QCD [27, 28]. 
In the present problem, Dalitz plots for — > if~7r+7r+ 



involve energies in the range 0.4 < Mf^ < 3 GeV, where 
the K is taken as particle 2, adopting PDG[29] conven- 
tion. The low-energy end of this range is directly within 
the scope of ChPT and it is in this region where discrep- 
ancies between LASS[1] and FOCUS [9] data are more 
pronounced. Evcin outside this range, the chiral frame- 
work still proves to be very useful and provides rigorous 
guidance for possible extensions. 

The case of the k is paradigmatic. This state was found 
in a number of experiments [3, 8-10, 30] and theoreti- 
cal models [3 1-34], with a complex mass = [(0.672 =f 
0.040)-i (0.550-0.034)] GeV[29], which can be consid- 
ered as low in the hadronic scale. For this very rea- 
son, it definitely cannot be accommodated as a funda- 
mental resonance in a chiral lagrangian[28]. In ChPT, 
resonances can only be introduced as next-to-leading or- 
der, which decay by means of explicit couplings to pseu- 
doscalar mesons. In this framework, a resonance should 
become a stable particle when its coupling to the Kw 
system is turned off. 

The K does not fall in this category. Its relationship 
with ChPT is a more subtle one and has been clarified 
about a decade ago[32], with the realization that unita- 
rization of low-energy chiral results gives rise to ampli- 
tudes which have poles in the complex s-plane. The k 
corresponds to the state with the lowest energy and it 
is present even if one considers only leading order con- 
tact interactions (see, for instance, our fig. 4). This idea 
proved to be very fruitful and motivated both reanalysis 
of low-energy experimental Kn data[33, 34] and the inter- 
pretation of the K as a tetraquark or two-meson state[35]. 
The case of the a- meson, which occurs in SU{2), is to- 
tally similar ( see e.g. [36-38]). There, however, the avail- 
ability of precision data allows the pole to be extracted 
directly from the tttt amplitude, in a model independent 
way [3 9]. 

Our presentation is divided as follows. In sect. II we 
review, for the sake of completeness, chiral results for 
the Kw amplitude and its unitarization. The conceptual 
framework for the three-body unitarization is introduced 
in sect. Ill and, in sect. IV, its perturbative expansion is 
used to derive predictions for the £>+ — )• K~'jt~^w~^ am- 
plitude. Results and conclusions arc given in sect.V. The 
manuscript also includes several appendices dealing with 
technical matters. 



II. Kir AMPLITUDE 
A. Interaction kernel 

The reaction WaiPa) Kb{pb) ndPc) Kd{pd) is de- 
scribed in terms of the usual Mandclstam variables s, t, u, 
constrained by the condition s + t + u — Pa+Pl + 
Pc Pd ■ results can be expressed in 

terms of the momentum q, with = sp^/A, p = 
VI - 2 {Ml + MD/s + {Ml-Mlf/s\ 
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Chiral perturbation theory determines the tree amph- 
tudes T/, with isospin /, as the sum of a 0{q^) con- 
tact term[27], supplemented by 0{q'^) scalar and vector 
resonances [28], together with inelastic contributions [33]. 
The ampUtude for each channel, indicated in fig.l, is then 
written as 



T, = T'^ 



(1) 



in eq.(2). The neglect of these terms allows the Ktt am- 
pUtude to be unitarized in a very compact form, as dis- 
cussed in the sequence. The 1 = 3/2 amplitude is re- 
pulsive, contributes little to FSIs and is also neglected in 
the sequence. Projecting out the S'-wave, one finds the 
leading kernel of the dynamical equation, given by 



K,=K,c 



N I s 
.1 — m? 



(5) 



FIG. 1. K-K amplitude: the contact diagram is the leading 
one at low-energies, whereas blobs are corrections due to other 
intermediate states; 5* and V correspond to scalar and vector 
resonances and I, to inelastic channels; the summation signs 
indicate the possibility of more than one intermediate state 
of each kind. 

In this exploratory work, we are interested in deter- 
mining dominant structures in the scalar sector and, in 
the sequence, we neglect vector resonances and inelas- 
ticities. As shown in fig.l, scalar resonances contribute 
to s, t and u channels. Each of the corresponding am- 
plitudes has the form T^{x) — aj,(x)/(a; — m^) , where 
x = s, M , is a ©(q^) polynomial and ttIx ^ a large 
resonance mass. In the physical region, the variables t 
and u are negative, whereas s is positive. This means 
that r^(s) can become arbitrarily large, whereas T^(£) 
and f^(u) remain always finite. Close to threshold, on 
the other hand, all the T^{x) are smaller, by construc- 
tion, than the contact term T"^ . Therefore, in the entire 
physical region, the combination \T^(t) -h T^(u)\ can be 
neglected in comparison with [r'=-t-T'5(s)] . We thus iden- 
tify the leading tree Kix amplitudes with the functions 



1 



F2 

3 
2F4 



-Mi)] 



N 

E 



-Mi) 



3/2 



1 

2F^ 



[s-(M2 + M|.)] 



(2) 



(3) 



(4) 



where F is the meson decay constant, N is the num- 
ber of scalar resonances of the = 0+ considered, with 
masses rrii and coupled to two mesons by the constants 



and 



c^, estimated in ref.[28] 



Their values are such 
that Q Cm and {cdCm/rnf) ~ 10"'^, explaining why 
f^{t) and f^iu) are much smaller than the terms kept 



/C,= [5s/8-(Af2 + M|.)/4 
-3(Mj 



(6) 



This kernel is real and hence suited for describing elas- 
tic processes only. The inclusion of inelasticities, due 
mostly to intermediate states containing rji and rjg,, can 
be performed by means of well known coupled-channels 
techniques [33]. We remain within the elastic approxima- 
tion and derive T, the elastic / = 1/2, S'-wave ttK scat- 
tering amplitude, by means of the Bethe-Salpeter equa- 
tion (BSE), written schematically as 



(2^ 



1 



[{e+p/2)^-M^+ie] 



[(£-p/2)2-Mi+z6] 



(7) 



with — s. Our approximations ensure that both /C and 
T do not depend on £ and the BSE acquires the very 
simple form 



T=[l-Tfl] K., 
where the two-meson propagator is 



(8) 



dH 1 

(2^ [{£+p/2y-M^+ie] 
1 



[(£-p/2)2-Mi+z, 



(9) 



This integral is ultraviolet divergent, but this unwanted 
behaviour can be cured by means of a subtraction. Fol- 
lowing ref.[27], we write the regular part of as 



n{s) = n{s) - n{o) 



(10) 



where the divergent part of Q{s) is contained in i7(0). 
Regularization amounts to the replacement 51(0) — ?> C, 
where C is an unknown finite constant which has to be 
fixed by experimental input. The function f2 can be eval- 
uated analytically, and explicit results are given in ap- 
pendix A . After regularization, eq.(8) becomes 



T=[l-T{C + n)] IC, 
and its solution reads 

T = IC/D , 

D^l + {C + n)IC. 



(11) 

(12) 
(13) 
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The Kit amplitude is thus determined by a rather 
simple algebraic equation, involving just the two-meson 
propagator and the kernel. The latter, in turn, includes 
only contact interactions and s-channel resonances. In 
the sequence we will use the result 



T {C + n) = i- l/D 



(14) 



derived from eq.(12). The function is real below the 
threshold at s = (M^+Mk)^ and acquires an imaginary 
component above it. In the physical region, below the 
first inelastic threshold, the kernel is real, whereas f2 — 
Clu + iClj, Qj — — p/IGtt, and the unitary amplitude can 
be represented by means of a real phase shift S [40] as 



T = 



167r 



sin (5 f 



tan (5 — — 



VLiJC 



l + (C+r!i?)/C 



.(15) 



For energies above inelastic thresholds, the kernel K, ac- 
quires imaginary components owing to processes included 
in the blob / of fig.l and the amplitude is damped. 



B. Alternative representations 

As a side comment, we note that the representation 
of T, in terms of Breit-Wigner structures, is not suited 
to this problem. They apply just to a single isolated 
resonance, since eq.(5) would read K. — — ai/(s — m\) 
and eq.(12) could be written as 



T = 



s — [m\ + {C + rifl) ai] -I- i [pai/lGTr] 



(16) 



In this case, the terms [m\ + {C -\-Q.p:) ai] and [pai/lQn] 
could be identified as a running mass and a width respec- 
tively. On the other hand, in the more realistic situation 
where both the contact interaction and other resonances 
are present, the coupled structure of the problem shows 
up. For instance, in the case of a contact term supple- 
mented by two resonances, the amplitude (12) would read 



(s — m^) (s — TO2) A^c~ (s — wi) Q^i ~ {s — m\) a2 
l + {C+nH + ini)[K.c- ai/{s-ml)-a2/{s-ml)] 

^ 7 2T7 2T • (-'-''') 

(S — mf) [3 — 1712) 

The coupling of various types of interaction gives rise to a 
complicated structure which cannot be written naturally 
as either products or sums of individual Breit-Wigner 
expressions for each resonance. 



C. Data and poles - schematic features 

The main qualitative features of the FSIs can be un- 
derstood by means of a Ktt amplitude as given just by a 
chiral contact term, supplemented by a single resonance. 
When more resonances are included, these features are 
preserved, but the amount of algebraic work increases 



considerably. We work with the simplest version, which 
corresponds to 



T = 



1 + (C + 17)/C ' 

/C = /Cc — oil/ {s~m\) 



(18) 



This amplitude depends on two sets of 3 parame- 
ters, namely [F,Mt^,Mk] and [mi,c^,c^], besides the 
subtraction constant C. The first set is determined 
by chiral perturbation theory and we adopt F = 
y/F^FK = 0.102722 GeV[33], M+ = 0.1396 GeV, 
= 0.4937 GeV. The other one is obtained by fit- 
ting LASS data[l] in the elastic region and we find 
C = -1 X 10-3, TOi = 1.330 GeV, ^ = 0.0352 GeV 
and c,\j — 0.001027 GeV, which yield the curve shown in 
fig. 2. Results for the coupling constants are close to one 
of the sets given in ref.[33] [cd,c„ = 0.030, 0.043]Gey 
and roughly consistent with those derived from the de- 
cay ao — rjTT in ref.[28], namely Cd = 0.032 GeV and 
= 0.042 GeV. 




Vs (GeV) 

FIG. 2. S-wave isospin 1/2 phase shift: LASS data[l] and our 
fit. 



Feeding the parameters from the fit into eq.(18) and 
extending it to the second Riemann sheet, one finds 
two coupled poles in the complex s-plane, given by the 
condition D{s) = 0. They are located at s — di, 
with = 0.51426725-^0.51423116 GeV^ and 6*1 = 
2.16256534-^0.24130498 GeV^, and identified respec- 
tively with the K and the ifg(1430). The latter is com- 
patible with the PDG[29] value, 5(1430) = [(2.01 ±0.15)- 
i (0.38 ± 0.13)] GeV^. In order to produce a feeling for 
the strength of the coupling between these resonances, 
we keep just the contact term in eq.(18), by setting 
ai = 0, finding 9^ = 0.45505779 - i 0.51167711 GeV^ 
and 9i = ml. The K-pole is thus rather stable. The be- 
haviour of the function D{s) along the complex s-plane 
is shown in fig. 3 and it is possible to see two wells around 
the K and Kq{14:30) poles. The plot in fig. 4 corresponds 
to the case ai = , in which just the k exists, dynami- 
cally generated. 

The fact that the function D{s) has two poles in 
the second Riemann sheet, allows the elastic amplitude. 
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FIG. 3. Two poles in \D{s)\ in the complex s plane. 




FIG. 4. A single pole in \D{s)\ in the complex s plane. 



eq.(12), to be represented as 



71 



(19) 



with /? 0.2200 , -0.1849 + ^0.6378 GeV^, 71 = 

0.2247 + i 0.1260 GeV^. In fig.5 we compare it with the 
exact amplitude and learn that this simple representation 
is reasonable in the range covered by the Dalitz plot. We 
adopt it in the present exploratory work and leave a more 
complete treatment for a future investigation. 



where W is the weak vertex and the functions a incorpo- 
rate hadronic FSIs. 

In figs. 6 and 7 we show the structure of A in terms 
of the rescattering series, which depends on T , the Kt: 
amplitude obtained in the preceding section. The dia- 
gram involving just W in fig. 6 describes the possibility 
that the mesons produced in the decay reach the detector 
without interacting. As the Fermi coupling constant Gp 
entering W is very small, the amplitude A(mf2,TO|3) is 
linear in this parameter, indicating that each dynamical 
component of the weak vertex in fig. 6 evolves indepen- 
dently by means of FSIs. In other words, if the primary 
vertex can be written as a sum W — of dynam- 

ical contributions, one automatically has A = v4(i) . 

In principle, tt'^tt'^ interactions could contribute to 
FSIs, but this subsystem has isospin 2 and this channel 
can be safely neglected, because phase shifts are small. 
Three-body unitarity is then dominated by the series rep- 
resented in fig. 7, where the first diagram corresponds 
to the spectator-pion approximation, often found in the 
literature. Contributions from three-body irreducible 
rescattering process are terms of the series which con- 
tain successive interactions among different pairs. The 
second diagram is thus the lowest order contribution of 
this kind. 



FIG. 6. Diagrammatic representation of the heavy meson 
decay process into Knn, starting from the partonic amplitude 
(red) and adding hadronic multiple scattering in the ladder 
approximation. 
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FIG. 5. Real (full lines) and imaginary (dashed lines) com- 
ponents of the Kn amplitude from eq.(18)(blue) and the pole 
representation (PR), eq.(19) (red). 



III. THREE-BODY UNITARITY 

The amplitude K^tt^tt^, denoted by 

A{ra\2, '7123), is symmetric under the exchange of the fi- 
nal pious and written as 



A{rn\ 



-W{m\2,ml^) + aiml^) + a{7nl^) ,(20) 



FIG. 7. Rescattering series implementing three-body unitary. 

In this work we concentrate on the perturbative struc- 
ture of the series describing strong final state interac- 
tions. In order to perform such a study, one needs in- 
formation concerning the production of the Xtttt system 
in the primary weak decay. A suitable conceptual point 
of departure is the quark diagram approach described by 
Chao[41], based on 6 independent topologies. In spite 
of their symbolic character and the absence of interac- 
tions mediated by gluons, they implement properly the 
CKM quark mixing. In the specific case of non-leptonic 
D decays, leading contributions are incorporated into the 
hadronic effective lagrangian given by Bauer, Stech and 
Wirbel[42] as 

G F 

icff = {fli {u d!)H {s' c)h + 02 (s' d')H {u c)h}(?1) 

where {qqjH represents hadronic expectation values of 
{V—A) charged weak currents and, the QCD factors al 
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and a2 are related to the Wilson coefRcients with ratio 
roughly given by 02/01 ^ —0.4. A detailed study of 
the process K^tt^tt^ based on this lagrangian 

has been performed by Boito and Escribano[25], consid- 
ering two-body FSIs only. These authors assessed the 
relative importance of contributions proportional to ei- 
ther fli or a2 to observables and the latter was found 
to rather visible, although systematically smaller than 
the former. In this work we focus on the strong evo- 
lution of leading contributions and neglect, for the mo- 
ment, the color suppressed term proportional to a2 . The 
light quark sector in the term proportional to ai involves 
just {ud)H 1 which is minimally realized by the matrix 
elements ( 7r+|A+|0 ) and ( 7r''7r+|y+|0 ). Concern- 
ing the factor {sc)h, we allow the final strange quark 
to be carried by either a kaon or a strange scalar res- 
onance, and arrive at the set of topologies indicated 
schematically in fig. 8. The strengths of these vertices 
are respectively Wa, Wb and Wc, assumed provisionally 
as constants. In the case of Wa, one h as t he isospin 
substructure: K~{k) n+{q) TT+{q') ^/2/3 Wa and 
K°ik) 7T°{q) TT+iq')^-^ Wa. 
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FIG. 8. Topologies for the weak vertex: the dotted line is a 
scalar resonance and the wavy line is a , which is con- 
tracted to a point in calculation; in diagram c, one of the 
pions is neutral. 

The construction of the amplitudes 0(^.^2) discussed 
in the sequence and formulated in terms of the rather 
general Kn amplitude given by eq.(12). Hence it does not 
depend on values adopted for parameters, the number of 
explicit resonances considered and possible couplings to 
inelastic channels. The function a(m^2) is written as 



a(m^2) = 



^=1 



(22) 



where N is the number of Kir amplitudes intervening in a 
particular diagram of fig. 7. The sum in (22) can be per- 
formed by means of a Faddeev-like decomposition of the 



amplitude A{r 



The Faddeev components are 



identified with a{mi2) and a{m2^), which are the non- 
perturbative solution of a scattering equation. 

Following the model proposed in ref. [16], we write the 
full decay amplitude of fig. 7 as: 



-4(mi2,m23) = W 



d*qd\'T3^3{k,k';q,q') 
(27r)8 (q^-M^+ie) 



(g'2 - Ml -I- i e) [{P -q' - g)2 - Mj^ + i e] 



(23) 



where P, k and k' are respectively the momenta of the 
D and of the pions produced in its decay. The matrix el- 
ement of the 3 — 3 transition matrix is Tj,^^{k, k'; q, q'). 



In order to simplify the description, we use the point-like 
weak vertex. 

The weak vertex for the decay of the D meson into 
the Ktttt channel is convoluted with the 3 — > 3 off-shell 
transition matrix, which takes into account the three- 
meson interacting final state, as shown in fig. 7, including 
the three-body connected ladder series, where the 2 — 2 
scattering process is summed up in the Kt: transition 
matrix. 

The 3 — >■ 3 transition matrix is obtained from the fol- 
lowing assumptions: i) the Kirn Bethe-Salpeter equation 
is solved in the ladder approximation, and ii) the Ktt 
transition matrix is effective in the S-wave states. The 
full 3—5-3 ladder scattering series is summed up when 
the integral equations for the Faddeev decomposition of 
the scattering matrix are solved. 

The matrix elements of our Kn amplitude depend just 
on the Mandelstam variable s of the two-body subsystem, 
as given by eq.(8) allowing the Faddeev components of 
the decay amplitude to be written in a factorized form: 



a(m?2) = Timj^} Cfe) 



(24) 



The function a(mi2) thus carries the full effect of the 
final state interaction through the two-meson amplitude 
T multiplied by a spectator amplitude ^, which contains 
the full three-body rescattering contributions, given by 
the sum of all processes in the ladder approximation for 
the multiple scattering series. 

The re-summation of the scattering series in the re- 
duced amplitude ^(fc) can be done by an integral equation 
corresponding to fig. 9, given by 



Cik) = ^,{k)~t 



d'^q T[{P-qf] 
(27r)4 (g2 - M2 + ie) 



[{P-k^ 

where the driving term is 



-g)2-M|. + ie]' 



~iW 



d'^q 



1 



(27r)4 (g2 _ M2 + le) 
1 



'~ [{P~k~q)^ ~ M]^ + ieY ^^^^ 
and depends just on a one-loop integral. This term, mul- 
tiplied by T{ 

1^12), gives rise to the first diagram in fig. 7, 
whereas the second term in eq.(25), multiplied by T{rn\2) 
represents the sum of all the remaining three-body con- 
nected processes shown in that figure. The lowest order 
connected contribution is the second diagram in fig. 9. 

The contribution to the three-body rescattering pro- 
cess given by eq.(25) is built by mixing interactions from 
the two possible Kn pairs. The resulting reduced am- 
plitude is a function of just the momentum of the final 
spectator pion. The model separates the decay amplitude 
into two parts: one corresponding to a smooth function 
of the momenta of the pions in W and another given by 
^(/c) times the pair amplitude, which is a fully three-body 
interacting term modulated by the Kn amplitude. 
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FIG. 9. Diagrammatic representation of the integral equation 
for the three-body function T{mi2)£,{k3) (left). The driving 
term contains the partonic amplitude from the weak vertex 
convoluted with the two-body scattering amplitude (right, 
first graph). 



So far, we did not consider isospin degrees of free- 
dom. The kernel of the integral equation (25) involves 
the change between the pions corresponding to the final 
isospin channel of the pair Kir and gives rise to isospin 
factors discussed in appendix B. The re-coupling coeffi- 
cient given by eq.(B3) appears weighting the kernel. Tak- 
ing into account also the isospin weight for the driving 
term, we find 



(g2 _ M2 



[{P-k-qf 



hp 



(26) 



As the main purpose of this work is to investigate the 
effect of the three-body unitarity on the decay amplitude, 
we analyze in the following section the perturbative con- 
tributions to the FSI at one- and two-loop approxima- 
tions. We choose to exemplify the series expansion of 
eq.(26) for the weak vertices a and c of fig. 8 and present 
the case of vertex b when discussing the perturbative cal- 
culation. In the case of vertex a, we find 



Aaiml 



2' ' 



^23) 



Wa\l + T{r 



'12 



42) 



2 ^ . 2 ^ 
+ 36^2)+ (3 

+ (1 o 3) . 



6("ti2) 



(27) 



where the argument of the function ^ is written in terms 
of the invariant mass squared of the Ktt subsystem, i.e., 
= (P—Pa)^, instead of the individual momenta. The 

factor comes from the isospin projection of the Ktt 

pair in the weak vertices to / = 1/2. The perturbative 
n-loop amplitude is constructed recursively as: 

A n{p-qr] 



X 



(2^)4 (g2 _ Af 2 + ie) 
Cn-l[iP-q)'] 



[{P-k- qf 



Ml 



(28) 



For later convenience we introduce the function 

2), 

2 \ rrt / „ _ 2 \ / 2 



Anl^io)' defined as 



A„(mi2) = T(toi2)'?»("1i2) 



(29) 



which is useful within our approximation of disregarding 
the momentum structure of the weak vertex, and eq.(27) 



becomes 



Aa(mi2,mi3) 



Wa{l 



I ) -^3(™?2)' 



Ai(m?2) 



A2(mi2)+ 



+ (1 O 3). 



(30) 



The three-body re-scattering series starting from the 
weak vertex b has to be treated properly in order to avoid 
double counting in the scattering series in the two-meson 
channel, as the scalar resonance is dressed by the Ktt 
interaction (c.f. fig. 12). In the case of vertex c the 
scattering series simplifies as the ttq produced directly 
from the W decay is not present in the final state and it 
is written as: 



Ac{ml2,ml 



V2 



Ai(m^2) 



A2(m^2) 



A3(m?2)--- +(1^3). (31) 



IV. PERTURBATIVE PROCESSES 



In this section, the first two terms of the function 
0(771^2) given by eq.(22) are evaluated covariantly and dif- 
ferent contributions are classified according to the type of 
initial weak vertex. Diagrams involve two kinds of loops, 
containing either two or three meson propagators. The 
former require regularization and are treated as in the 
construction of the Kn amplitude presented in sect. II. 
The latter are triangle integrals, written as 



P. 



KB 



(W12) 



1 



(27r)4 [(pi2-fc)2-Af2 + ,,] 

1 

[[k^-Ml+ie] {p3 + k)^-e] 



(32) 



where 6 = 9^ — iOi, is the position of the pole in the 
complex s- plane, with 0^ and 9i constant positive quan- 
tities. This integral is similar to those occurring in usual 
calculations of form factors, but not identical, since the 
invariant masses along the dotted lines in fig. 10 can be 
smaller than either to|j or m^2- It is thus mathematically 
more akin to integrals needed to describe form factors of 
unstable particles, such as the A, p or K* . We write 



/^A'e = zn^ife/(47r)^ 



(33) 



and the evaluation of the functions 11 is discussed in ap- 
pendix D. 



A. Contributions proportional to Wa 

Processes involving the weak vertex Wa, defined in 
fig. 8, are indicated in fig. 11. The W'^ is shown explicitly 
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FIG. 10. Triangle diagram representing the reaction D{P) 
[■kK]{pi2) vr(p3), with intermediate states 7r(pi2— fc), 0{pii+k), 
K{k), associated with the integral given in eq.(32); invariant 
masses along the dotted lines can be smaller than those of 
external lines. 



4^ 



FIG. 11. Diagrams involving the weak vertex Wa', the wavy 
hne is a , always plugged to a tt'*'; the tt produced to- 
gether with the K in the opposite side can be either positive 
or neutral. 



just for the sake of clarifying the various topologies, and 
is taken as a point in calculations. 

We start with the perturbative expansion given in 
eq.(30), and we evaluate the one and two- loop terms 
as follows. Using eqs.(25) and (9), one gets Ci("^i2) = 
-Wail(mi2) and Ai(mf2) = - r(m^2) ^^("^?2)■ As is 
divergent, a subtraction is needed and we assume the sub- 
traction constant to be the same as in the Kn scattering 
amplitude. Using (14), one has 

= -l + X[{ml,), (34) 



(35) 



Our assumption allows further simplifications in the 
treatment of the rescattering series, but it is still a free- 
dom within our framework that we will not explore fur- 
ther, in order to minimize the number of free parameters 
in our first investigation of three-body rescattering ef- 
fects. 

The two-loop contribution comes from the recursive 
formula (28) and eq. (29), and reads 



1 



(27r)4 [(pi2-fc)2-M2 + ze] 



1 



(36) 

and, again, the divergent function ft shows up. Subtract- 
ing, we have 

A2(mi2)=«r(TOi2) 



1 



(27r)4 [(pi2-fc)2-M2 + ze] 



X Jt(P3 + fc)^]^ {C + n[ip, + kr]}. (37) 



[fc2 



Using eq.(14) and subtracting once more, we get 

A2(m?2) = T(to?2) [C + n{ml,)] + y,{ml,) 

= -Xi{ml2) + X'2{ml2) , (38) 



d^fc 1 

(2^ [{pi2-ky-M^ + ie] 

^ (39) 



Ml-fze] [D(p3 + fc)2] ■ 



Repeating this procedure for higher order loop contri- 
butions, we find the structure 



XNiml^) = -AAr_i(mi2) + A^(toi2) 



(40) 



which can be checked explicitly for A3 and A4 , using the 
expressions shown in appendix C. Denoting by S the 
sum of terms within the square bracket in eq.(30), we 
have 



Af-l 



N=l 



til 



N=l 



N-1 



A', 



N 



(41) 



This allows eq.(30) to expressed as 



Aa{m\2,ml^) 



A'i("^i2) + 3 Ki^nl^) 



+ (0'a3(™?2) + 



(1 o 3).(42) 



This result is interesting and indicates the importance of 
treating both the Kir amplitude and the FSIs in a single 
coherent dynamical framework. The diagrams of fig. 11 
contribute to the first two terms of this series. 

The integral (39) can be performed by a number of dif- 
ferent techniques, provided one recalls that 1/D{s) con- 
tains two poles, associated with the n and the if* (1430). 
One possibility would be to perform a Cauchy integration 
over fco, supplemented by a numerical three-dimensional 
integration. An alternative, adopted here, is to rely on 
usual Feynman techniques. With this purpose in mind, 
we employ the expression for the Kir amplitude written 
in terms of its poles, given in eq.(19) and, using (32), 
obtain 

A2(m22)-T(m?2){-/3 [C+^iml^)] 

+ 7«n,/fe J(167r2) + 710^x61 /(levr^)} (43) 

The problem then reduces to the evaluation of triangle 
integrals. 
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FIG. 13. Diagrams involving Wc', one of the pions in the weak 
vertex is neutral. 



FIG. 12. Diagrams involving the weak vertex Wi,; the wavy 
line is a , always plugged to a vr'*' and the dotted hue is a 
scalar resonance, which has a width given by the substructure 
R described at the bottom hne. 



B. 



Contributions proportional to Wb 



series as in eq.(30). Using result (41) one finds 



-1 + A'i(m22) + | A'2(m?2) 



+ (|)'a3("^?2) 



(1^3). (49) 



We follow here the treatment of the production ampli- 
tude for a scalar resonance given in ref.[43]. The basic 
idea is that processes, shown at the bottom line of fig. 12, 
must always be considered together, since the tree contri- 
bution in isolation contains a pole in the physical region. 
The corresponding amplitude W, involving a resonance 
with mass rrii, reads 



[1-T{ml^)n{ml^)] . (44) 



Regularizing the function D, and using eq.(14) one finds 



K2 - mf] D{ml^) 



(45) 



This function is finite because, by construction, D(rnf) = 
0. 

The evaluation of the contributions to Ab is straight- 
forward and yields 



Ab ~ abi + ab2 + ■ ■ ■ 



(46) 
(47) 



(27r)4 [(pi2-fc)2-M2+ze] 

(48) 



where the recursive formula (28) is used starting with 
f 1(771^2) = fl6i(™i2)- The function af,2 can be reduced to 
a sum of triangle integrals and is evaluated using eq.(32). 



C. Contributions proportional to Wc 

Processes given in fig. 13 give rise to eq.(31), which con- 
tains the same 



V. RESULTS, CONCLUSIONS AND SUMMARY 

In this work we studied the role of final state inter- 
actions in i4r^7r+7r+, treating two- and three- 
body interactions in a consistent way. Our underlying 
Ktt amplitude is derived from chiral perturbation the- 
ory, supplemented by unitarization and tuned to elastic 
LASS data[l]. Two poles, associated with the k and the 
A'* (1430), were then determined and employed in a rep- 
resentation of the amplitude. A number of simplifying 
assumptions were made in the Kn amplitude, for the sake 
of minimizing technical problems. Among them, we men- 
tion the absence of both isospin 3/2 and P waves, as well 
as couplings to vector mesons and to inelastic channels. 
The means for correcting these shortcomings are avail- 
able in the literature and will be considered in future 
extensions of this work. Our treatment of three-body 
unitarity departs from a Faddeev-like integral equation, 
which is subsequently expanded perturbatively. Terms in 
the corresponding series contain a recursive component, 
which allows a re-summation of the whole series, when 
the divergence is subtracted using the same criteria as in 
the Ktt scattering amplitude. 

Model independent analyses [8- 10] of the s-wave if"7r+ 
channel in the decay D"*" — >■ i^T ~7r+7r+ are rather welcome 
in theoretical studies, because they are expressed in terms 
of amplitudes which are linear in the Fermi constant. 
This means that it is meaningful to study independently 
the strong evolution of each dynamical contribution to 
the primary weak vertex. In this work we have consid- 
ered just three simple topologies of the color-allowed type 
to the primary vertex, which yield classes of decay am- 
plitudes denoted by Aa , Ab , and Ac eqs.(42, 47,48, 49). 
The first two begin at tree level and give rise to quasi 
two-body FSIs, as in ref.[25]. The amplitude Ac, on the 
other hand, arises only when proper three-body FSIs are 
present. 

Predictions for their modulus are given in figs. 14. It 
is important to note that the order N of partial contri- 
butions indicates the number of times the denominator 
1/D oi eq.(14) intervenes in a given function. Results 
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for Aa and Ac, in particular, are based on an infinite re- 
summation of terms, eq.(41), and hence N is not simply 
related with perturbative counting. Understood in this 
sense, second order terms tend to be smaller that leading 
ones. However, in the case of Aa and Ac, FSI corrections 
show up clearly around the bare resonance mass, where 
the function 1/D vanishes along the real axis. In the case 
of Ab, one notices a cancellation close to threshold. As 
far as comparison with FOCUS data[9] are concerned, we 
note that Ac has a dip at the correct position, whereas 
compatibility with the direct production of a resonance 
at the weak vertex is very difficult. 
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m,,{GeV) 



2.0 

1.5 

I 1.0 
■o 
o 
E 

0.5 
0.0 
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0.8 



1.2 1.6 
m,, (GeV) 



FIG. 14. up: behaviour of \Aa\ and \Ac\; down: behaviour of 
\Ab\; first and second order partial contributions are indicated 
by dashed and dotted lines. 

Predictions for the phase are displayed in fig. 15. Par- 
cial contribution from X'l in eq.(42) and eq.(49) and from 
ttbi in eq.(47) fall exactly over the elastic Kit phase. The 
oscillation in Ai, at low-energies can be ascribed to the 
lack of precision in ab2 ■ We have checked that it is due to 
an incomplete cancellation between a precise tree term 
and the less precise triangle contribution in processes 
shown in fig. 12. Finally, we see that a curve for Ac, 
shifted by -148", describes well FOCUS data[9] up to 
the region of the peak. 



As the study of both the modulus and the phase seems 
to favor the weak vertex Wc of fig. 8, it is worth explor- 
ing its structure. In fig. 16, we show the phases of the 



factors [-1 + Ai(m^2)] and [-1 -I- X'liml^) + | A2(m^2)] 
in eq.(49). The gap between them is about 10" and the 
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FIG. 15. Phases from Aa , At , and Ac, continuous line; par- 
tial contributions from A'l in eqs.(42) and (49) and from an in 
eq.(47) coincide with the elastic Ktt phase, dashed lines; con- 
tributions from X'2 and ab2 are given, respectively, by dotted 
and dashed-dotted lines; the curve coinciding partially with 
FOCUS data[9] (triangle) is Ac, shifted by -148°. 



first term already corresponds to a good approximation 
to low-energy FOCUS data. However, as pointed out 
above, the factor —1 in these terms comes from an in- 
finite re-summation and hence the smallness of this gap 
does not indicate that it is enough to consider just the 
first perturbative contribution. Using eqs.(49), (35) and 
(12), one has 



Ac{ml 



27 "'23 



+ (1 o 3) 



-1 



1 



i + ic + n)ic 



(50) 



where /C is the Kt: kernel of eq.(5). Hence the dip in 
experimental results for the magnitude of A indicates the 
energy at which IC = 0. In our model, this corresponds 
to the point where both the real and imaginary parts of 
the scattering amplitude cross the real axis in fig. 5. 
In this framework, eq.(50), rewritten as 



A r 2 2 N V2 ... 
Ac{m-^2,m23)~— Wc 



1 + {C + Q)K. 



+ (lo3) ,(51) 



could provide a much more convenient structure to be 
used in experimental analyses, since $1 is a well known 
elementary function. 

As stressed before, this work concentrates at tracking 
leading effects, and a number of important issues were 
left untouched. Among them, we mention the absence 
of P- waves, isospin 3/2 terms and inelastic contributions 
to the basic Kn amplitude, as well as a more complete 
description of the primary weak vertex. These matters 
will be dealt with in other papers. 

In summary, our main conclusions read: 

1. proper three-body effects are important; 

2. corrections to the re-summed series are important at 
threshold and fade away rapidly at higher energies; 

3. a model based on a vector weak vertex can describe 
qualitative features of the modulus from FOCUS and 
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FIG. 16. Phases of the factors [-1 + Ai(m?2)] (dashed) 
and [-1 + Ai(m?2) + | AaCmia)] (continuous) in eq.(49), 
shifted by -148° compared with FOCUS [9] (triangle) and 
E791 [8] (circle) data, together with elastic K-k results from 
LASS [1] (diamond). 



as 



^{s) = n{s) - n{0) = - [L{s)-L{0)] /167r2 
= -L(s)/167r2 , 



where, 

.s<(M,-Mk)2 

Z = L + ^ In 



+ Ml - s + Va' 



{M^-MKf<s<{Ml+M]i) 



L = L- 



tan 



Ml + Ml - s 



• (Ml + Ml)<s<{M^ + MKf 



L = L- 



tan 



M2 + MI-S 



(A3) 



(A4) 



(A5) 



TT (A6) 



E791 data[8, 9] for the decay amplitude and agrees well 
with its phase in the elastic region. 
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.s>(M^ + MKf 

L = L- ^ In 
s 



+ in 



Va 



s-M^- Ml + ^/X 
2 Mk 



with 

L = l + 4±^ln 

M2 - Ml 



M^ 



Mk 



Ml - Ml 



In 



(A7) 



M^ 



Mk 



(A8) 



\^s' -2 s {Mi + Ml) + [Ml - MlY . (A9) 



Appendix A: two-meson propagator 



Appendix B: Isospin in FSIs 



FIG. 17. Bubble loop diagram. 



The basic integral is 



I-nK{s) 



(47r) 



(Al) 



since = s. In the Bethe-Salpeter equation, it is con- 
venient to use 



(A2) 



where L(s) is a finite function and A^o incorporates the 
ultraviolet divergence. The regular part of is defined 



In order to derive the isospin weighting factors pre- 
sented in the integral equation (26), one has to realize 
that the adopted model for the Ktt transition matrix has 
a separable form, which allows for the factorization of 
the decay amplitude as given by eq. (24). In addition, 
we consider only the dominant isospin 1/2 channel of the 
Kn system. The isospin 1/2 dependence in the transition 
matrix is made explicit by writing it in the form: 



ri/2 = J2 \^Kn = 1/2, I,) T {iKr. 



1/2, ^. 



(Bl) 



where the matrix element of T is given by eq. (19). 

The Kirn scattering amplitude with 1^ — 3/2 is built 
only with total isospin states It = 3/2, because the Ktt 
subsystem interacts only in the isospin doublet channel. 
The states with isospin 5/2 do not contribute to the final 
state interaction in our model. The above simplification 
implies that the spectator function ^, which also carries 
the bachelor pion momentum distribution, is given by: 



16/2 (fc.)) 



\It, Iz,Iir,lKTT' 



(B2) 
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where It = 3/2, 1^ = 3/2, I„ = I^, = \ and Iktv = 
W = 1/2. 

Following the diagrammatic representation of the inte- 
gral equation (26) shown in fig. 9, one reahzes that the 
kaon is exchanged between the two possible interacting 
Kit pairs. Guided by the diagram, one; verificis that the 
bachelor pion in the intermediate state interacts with the 
kaon and forms a new pair. This is a known property of 
the Faddccv decomposition, which equates only different 
Faddeev components of the 3— >-3 transition matrix. 

The kernel of the integral equation exchanges the role 
of the bachelor pions from the intermediate to the final 
state (see fig. 9). In addition, by considering the 1=1/2 
transition amplitude (Bl) and the isospin structure of 
the spectator function (B2), one finds that the weighting 
isospin factor multiplying the kernel of eq. (26) is given 
by the re-coupling coefficient 

-R3/2 = {It,Iz,I'!v,Ik7v'\It,Iz,I-!v' ,lK-!r) , (B3) 

which is written in terms of Wigncr 6-j symbol, 

R3/2 = i-lf f (Ik.) f {Ik.') { } ,(B4) 

with S = 2I^ + Ik + It and /(/) = ^21+1. The result 

is the weight factor i?3/2 = 2/3 in (26). 

In the particular case of eq. (26), the driving term for 
It = 3/2 is derived from a weak vertex symmetric under 
the exchange of the pions, such that: 

\Ws/2)=W{\lT,h,I.,lKn') + \lT,h,In',lK.)) . (B5) 

Both terms contribute to the driving term isospin weight- 
ing factor after projection onto the appropriate isospin 
state, which corresponds to the spectator function (B2). 
One gets that 

{It,Iz,I.,Ik.'\Ws/2) = 

= W{It, Iz, I-iT, Ikt7'\It, Iz, Iv, Iktt') + 

+W{It,Iz,I.,Ik.'\It,Iz,I.',Ik.) = Iw , (B6) 

which was computed by using the isospin re-coupling co- 
efficient (B4). 



Appendix C: Perturbative contributions 

Perturbative contributions to the FSl amplitude 
^("^12) — J2^=i «iv('7ii2) ; eq.(22), have the structure 

«Jv(™i2) = -ajv-i(mf2) + ai(mi2) Ajv(toi2) (CI) 
and, in the sequence, we list the A^(mf2) for A'' = 3,4. 

d-^k 1 



A3(Wi2)=-T(mi2) 



1 



(47r)4 {p^^-kf-M^+ie 

T{k+P3) 



d^k' 



(47r)4 [(/c+p3)-fc']2-M2+ie 
1 1 
^ k"^-Ml.+ie D[{k'+pi2-k)'^] ' 



(C2) 



X 



r d'^k' 1 

J (4^ [{k+p3)-k']^-M^+ie 

f d^k" 1 

J (4^ {[k'+{pi2-k)]-k"]y-M^+ie 
1 1 

X -77 



X 



fc"2-M|.-Fie D[{p3 + k-k' + k")^] 



Appendix D: Triangle integral 



(C3) 



The triangle integral defined by eq.(32) and repre- 
sented in fig. 10 is written as 



I. 



(47r)- 



(47r)^ 



da I 



db 



K0 



(Dl) 



D„Ke = {l-a)M^+a{l-b)Ml, + abeR 
— i[e+ab{6i — e)] — a{l — a){l — b) 
-a{l-a)bMl-aH{l-b)M^ . (D2) 

The double integral (Dl) can be evaluated numerically 
but, in the case 0i = e, problems of accuracy may arise. 
In order to understand the structure of its complex part, 
it is desirable to pursue the analytic path as long as pos- 
sible. We therefore resort to the SU (2) chiral limit and 
neglect M^, eliminating terms quadratic in b in Dsk-k 
and simplifying the algebra. Integration in b yields 



'-f 

Jo 



db 



G+i{9i-e) /I [F+G]2 + [(ej-e) + e/a]2 



G2 + (^j-e)2 \ 2 
. F{6i-e)-G{e/a) 



G{e/a)\ 



-1 x—y 

tan h cr TT 

l+xy 



F = Ml - (1-a) m 



2 

12 1 
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G={eR-Ml) - (1-a) {Ml-ml^) , 
(F+G) G+iOi-e) {Oi-e+e/a) 



x=- 



\F{ei-e)-G{e/a) 



y-- 



FG+jOi-e) (e/a) 
\Fi9i-e)-G{e/a)\ 



(D3) 



with [xy > —1] — )■ cr = , [xy < —1 and x > 0] ^ a = 
+1 , [xy < —1 and a; < 0] ^ a = — 1 . This result was 
used to tunc numerical ealeulations and one found that 
convergence for small values of e requires a large number 
of points in a Gaussian integration. 
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